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*ABSTRACT - pe 


The he eG 1S tracker system of optimal smoothing was investi- 
gated in a track-while-scan radar by using three dimensional simulation of 
various interceptor-target profiles. Random noise was put on all parameters 
which were not positively defined. This included the target's motion. 

The computer program used is designed for carrying on the analysis 
with other systems so that their performance can be compared with the 
=; & , 1 tracker. A number of representative profiles are listed 
which should provide a suitable basis for evaluation of any chosen 
alternate system. 

The parameters used in the analysis are similar to those which 
will be used in the Hughes Aircraft Co. developed AWG-9, Phoenix, 


Aircraft Missile Control System designed for the F-111B aircraft. 
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ee Introduction 

They = oe fe) sampled data system has been the most often 
proposed means of performing the task of optimum smoothing in a 
track-while-scan radar. This method is particularly attractive since 
it is very simple. However, in the light of the high performance intercept 
problems which may be encountered today, especially by a high perfor- 
mance aircraft, it was felt advisable to investigate the capabilities of 
such a smoothing system in the environment of air-to-air intercepts. 

In this analysis the third order system or ox P (S ; i 
tracker was considered by using a three dimensional simulation of various 
interceptor-target profiles. Random noise was placed on all the para- 
meters which were not positively defined. This included the motion of 
the target. 

The computer program used was deSigned for carrying on the analysis 
with other systems so that their performance could be compared with that 
of the Xx, & ; ae tracker. 

The parameters used in the analysis were similar to those which will 
be used in the Hughes Aircraft Company developed AWG-9, Phoenix, 
Aircraft Missile Control System. This system is currently being developed 


for use in the F-111B aircraft. 
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Zu Basic System Description 

The system considered in the analysis was a pulse doppler radar 
operating in a track-while-scan mode, TWS. That is, the radar was 
able to obtain simultaneously at discrete sampling intervals, called 
frame times, the target's range, range rate, and bearing, in azimuth and 
elevation, relative to an interceptor’s antenna. At the same time the 
system was also maintaining a track file on a number of other targets 
at locations spaced over its scan volume. 

The heart of a system such as this is the so-called "track-while- 


scan loop." Fig. 1 is a basic block diagram of this loop. 

Prior to entering the TWS loop the raw A to D converted parameters 
are taken from the output of the radar, and the interceptor's inertial 
reference platform and air data computer. This information is pre- 
processed to remove anyredundanciescaused by picking up the same 
target on two or more antenna sweeps, called bars. Also the north, 
east, down direction cosines are computed. TheSe are the direction 
cosines fora north, east, down axis coordinate system which is fixed 


to the interceptor. The axes are defined by having the north axis pointing 


toward true north. The east axis results from the cross product of the 
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down and north axis. The down axis points towards the earth's local 
vertical, and is positive downward. Thus an interceptor flying north 
would have the north axis along the line of sight, the east axis out the 
right wing and the down axis vertically downward. 

The data now enters the "track-while-$can loop." It is first 
correlated by having the individual observations checked against previous 
tracks on file in memory. This may result in several observations being 
linked with one track. These are then associated by having the observations 
paired by a logic routine with the track which it most closely matches. The 
next stage initiates new tracks or deletes old tracks, which no longer 
associate after a quality counter, which rates each track, degenerates to 
a point such that the track no longer warrants being maintained on file. 

The next section and the one which is investigated in this paper is 
that of smoothing and prediction. The ability of the system to,perform all 
previous and following tasks rests on how well the track is smoothed, i.e., 
the effect of noise removed, and how well it predicts the target's future 
motion. This section must be designed to respond quickly to target 
maneuvers while filtering out the false excursions of the target caused 


by noise so that future predictions will be as accurate as possible. 
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The gating section is next in the loop. It is this section which 
determines the possible limits of target motion between contacts. It allows 
for target maneuvers between scans plus approximately 30° noise. It further 
provides for expansion in the event of a missed contact. A maximum gate 
size is imposed which corresponds to a volume which would contain any 
possible target maneuver that could occur up to the time the track would 
normally be eliminated by the quality counter due to missed contacts. 

The individual track parameters are stored in the track file section 
to keep it updated. These along with the track locations are used to 
determine target priorities. The track parameters are also used to determine 
which observations correlate with the tracks during the correlation phase 
at the beginning of the loop. The above computations are performed on all 
tracks at the end of each frame time. 

Additional functions which result from the TWS loop computations 
are the antenna scan center calculations to keep the pattern centroid 
located in the center of the incoming raids, the steering information to the 
pilot to allow him to maneuver as necessary to keep the radar within its 
gimbal limits, prelaunch information to the missile so that its seeker head 


is oriented in the proper position at the beginning of its semi-active phase, 
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and midcourse maneuver corrections to the missile during flight. 

It is readily apparent that the system must be optimized to the highest 
possible degree in light of the high relative speeds encountered and the 
number of targets to be tracked. At the same time it is extremely important 
that computations remain relatively simple to maintain reliability and to 
keep computational time and storage space to a minimum. 

The ability of the system to meet these criteria lies primarily in the 
smoothing, prediction and gating sections. It is the former two which are 


considered in this paper. 


) b >’ vo 

b b , d be RD 
b ’ 
b 


wv 2se we ws ww 
od 
= 


d 


6 of 100 


aot tee a 4an> 
7 


7 6 
77 
ee see ot 





COST . 
’ % ‘ ¥ y 
’ ’ ’ ’ ’ y ¥ , »» 
y ’ 


% 
’ 
’ ’ 
’ 


‘ 
’% PO +vd)d ¥ ’ ,v} ¥» ¥ 


3: System Characteristics 

The values of the system performance parameters were derived 
primarily from the Hughes SYSTEM EVALUATION AND ANALYSIS REPORT, 
SEAR, /1 ni Where possible these parameters have been crossed checked 
with other company internal documents, or by conversations directly with 
company engineers. The following parameters and assumptions are those 
which were used in the simulation. 

It was assumed that the interceptor would have an on station Mach 
number of 0.7 at above 36,000 feet. The target's speed was believed to 
be between Mach 0.4 and 3.0 and its altitude could be between sea level 
and 80,000 feet. 

The general characteristics of the radar were those of a high PRF 
pulse doppler using linear frequency modulation. 

The detection range was 165 n.m. This was defined as the range 
where S/N was ore using the classic radar range equation. It was based 
on a five Square meter target. The radar's probability of detection curve 
is shown in Fig. 2. 

The effect of clutter on detection probabilities was not considered 
since a large amount of clutter rejection is an inherent advantage of a high 


PRF pulse doppler radar. This is particularly true over water. In addition 
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the system will also have clutter rejection circuitry. Due to the high PRF 
certain relative velocities may cause periodic loss of signal due to 
"“eclipising", i.e., target echo returning during pulse transmission. 
Multiple PRF rates will normally remove this possibility. Thus eclipising 
was not considered to effect detection probability. 

In the radar, range is measured by taking the difference between 
two filters in a doppler filter bank that are "rung" on successive FM phases. 
The spacing of the filters corresponds to five n.m. Since there are 30 
filters the range can be measured to 150 n.m. The signal was considered 
to originate uniformly in plus or minus one-half a filter and terminate 
normally distributed with a standard deviation of one-third a filter. 


This is shown in Fig. 3. 
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Hughes stated that a joint distribution describing these two distributions 
yielded a normal distribution of zero mean and standard deviation of 
1.37 n.m. It was felt that this value was too optimistic and indicated that 
the distributions were quite dependent, i.e., the covariance was large. 
The value used for the range standard deviation in this analysis was 
Zo) hn « 

Range rate is measured by frequency shift as determined by a doppler 
filter bank. The filters are spaced 150 cps apart. Thus at X band using 
the equation given by Skolnik /2/: 


V = 4 ~ 
. 103 


Where: f 


d doppler frequency shift. 


4 
il 


wavelength in cm., 3cm at X band. 


< 
li 


range rate in knots. 


The total number of filters is 768 giving the ability to measure range 
rates of up to approximately 3250 kts. 
It was determined by Hughes that the distribution of the range rate 


error was normal with zero mean and a standard deviation of two filters, 
B.f2 ktee 
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The scan pattern of the antenna in the TWS mode is either two bar 
+ 40 deg. or four bar + 20 deg. The beam separation is 1.32 deg. and the 
scan rate is 40 deg/sec. It was determined from various sources that the 
average antenna azimuth error could be considered normally distributed with 
zero mean and standard deviation of 0.745 deg. Similarly the elevation 
error had zero mean and standard deviation of 0.86 deg. These values 
include radome, potentiometer output, and pick off errors. The additional 
error in elevation is due primarily to the introduction of redundancies caused 
by picking up the same target on two Separate scans. 

The smoothing and prediction system which has been proposed by 
Hughes isthe @&X<, Q ; ‘ tracker type. The characteristics of this 


system are considered in the following section. 
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4, Smoothing System Discussion 

Since computational time and storage space available for smoothing 
and prediction is quite limited some of the more sophisticated smoothing 
systems would not at first appear to be suitable for use in the system 
described. The requirements of the system are reduced somewhat, however, 
by the fact that the target's capabilities are relatively predictable, but of 
course, entirely random. For these reasons the most attractive choice is 
an @—&, S tracker type of system. Almost any other system that would 
be acceptable would reduce to a scheme which would use a similar smoothing 
constant gain technique in any caSe. 

The S&S , & ; Ss tracker system is an extension of the xX ; ~ 
tracker to allow for the smoothing of acceleration. As given by Benedict 
and Bordner V3, the variance reduction ratio is optimized by having the 


velocity constant, iS , defined as : 
ee 
> 


Oa 


Similarly it was shown by Simpson /4/ that when acceleration was 
considered the optimum value of the acceleration smoothing constant, ¥ , 


became: 


ye ES - (x +8). 
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Thus if it is possible to define S&S the smoothing equations 
become unique. 

The proposed equations used for the smoothing and prediction 
are: 


A. Range: 
= - Ret o CR a~ Ren) 
ae Wes oT (Rr- Rew) 
. oo eae 
R Pad Ran + T Ry ae 


At 


Raw? ie * T Re 


Where: 


e 


ao and Tse -R ©, if an observation 
was received. 


eae and RR ys oe if an observation was missed. 


| =Frame Time = 2+ 0.1 Sec. 
sine Number of frame times since last value of range rate was 
received. 
The subscripts are: 
(\= Frame number. 
P= Predicted. 
© = Observed. 
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The value of “ _ for range smoothing as proposed by Hughes was 


0.45. Since the standard deviation of range rate was relatively small 


no smoothing of range rate was necesSary So that no © OG N 


were required in range smoothing. 


B. Direction Cosines: 


a) oC, Ny.) 

- ~A. . i ) + 2h he) 
= Ne. ) a (Nn Ne 

oe tamed he 


Where: 


_N\. oe aXe oy, if an observation was received. 


Ny DS >, if an observation was missed. 


| \. = general direction cosine being smoothed. 
For the North and East direction cosines, N, E, Hughes proposed 


the smoothing constant values of: 


eC 


S = 0.120 
ee 0.008. 
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For the down direction cosine, D, the proposed values were: 


SV = 0.35 


(© = 0.075 
O= 0.0035. 


A Z~ transform analysis of the & , e Tracker was performed by 
Sklansky Toys Basically the system is stable and for a non-maneuvering 
target there would be no steady state error for % between zero and one. 
When ™ is equal to one no smoothing would be performed and raw data 
would be the output. Ina high noise non-maneuvering target condition 
the smoothing constant would be small and the system response would 
be slow. For increasingly lower noise the smoothing constant would 
approach one in the limit. For a maneuvering target the smoothing 
response must be increased so that a larger value of &y is required. 

The obvious disadvantage of the >¥ , S Tracker is that there is 
not a single optimum value of the smoothing constant for all the situations 
that can be encountered. 

The advantage of the system again lies in its simplicity and the 
small number of computations and amount of storage space required. 
eS since © and G are defined uniquely when S< is determined the 


problem of fixing the constants is reduced somewhat. 
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Sig Simulation 


The CDC 1604 computer was used to implement the three dimensional 
simulation of a single target-interceptor problem. This program was named 
Three D. The primary aims of the simulation were to provide realistic, 
flexible and accurate target-interceptor geometries. The pseudo random 
noise generator, sub routine RNDEV, was used to provide gaussian noise 
on all parameters that were not positively defined. 

A list of symbols and a description of the use of the programs used 
are contained in Appendices A and B. 

The resulting program allowed the target to be positioned initially 
at any location in the northern half of a north, east, down coordinate 
system which was ground stabilized at the interceptor's initial location 
at problem time zero. The target's initial coordinates , maneuver para- 
meters, and maneuver times were then read in by data cards along with 
the interceptor's heading, altitude, and Mach number. The interceptor’s 


velocity was computed in the program by using the following relation: 


Vins = Sonic Velocity fps x Mach Number. 


The variation of the speed of sound with altitude is shown in Fig. 4. 
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Speed of Sound Versus Altitude 


No provision was included to maneuver the interceptor since one of 
the basic assumptions of the system was that the interceptor was continually 
tracking Several targets in widely spaced locations in the radar’s scan 
pattern. It was thus not possible for the pilot to place an individual target 
in the center of the scan. Also no random noisSe was uSed to vary the inter- 
ceptor's position since the same relative effect was gained by randomly 


varying the target's velocity components. 
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The problem time was taken to be 300 seconds. This allowed a target 
at Mach 1.5 above 36,000 feet to travel slightly over 70 n.m. in the course 
of the problem. Since the probability of acquiring and tracking a target 
beyond 100 n.m. was very low, 100 miles could be considered the range 
limit in setting up the problem profile. 

The target's velocity was also determined by the previously mentioned 
relations. No provision was made to accelerate the target, but the target 
would decelerate at the start of its dive. Its final speed was a function 
of altitude below 10,000 feet. The minimum speed was Mach 1.2, 1340 fps. 
at sea level. The relations used are shown in Fig. 5. The target was 
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onsen to make a 180 degree turn er a geen time and horizontal 
turn g load. It would also make a 60 degree dive and pull out at a pro- 
grammed time, final altitude, positive and negative g factor. 

For convenience, distances were read in, in nautical miles, velocities 
in Mach number, and headings in degrees. The program then converted the 
values to compatible units of ft., ft./sec., and radians. 

The motion of the target was computed by piece-wise—linear integration. 
Random accelerations having a standard deviation of one-half g for two second 
periods were imposed. The target's motion was broken down into velocity, 
VT, along its path of flight, heading angle, AT, and pitch angle, BT. Each 
component was then related to the distance traveled in each of the coordinate 
axes in the iteration time and combined with the previous position to give 
the new target location, XT, YT, ZT. 

Every two seconds of problem time the target's coordinates were taken, 
compared with the interceptor's coordinates XO, YO, ZO and the actual 
antenna azimuth, ATA, elevation, EPSI, range, R, and range rate, RDOT 
were computed. If the target had passed out of the limits of the radar beam, t 
65 degs. in azimuth and elevation the problem stopped. It should be noted 
that the antenna does not have a scan volume this large in the track while 
scan mode, but it does have the capability of positioning its scan anywhere 


inside these limits.. a i 
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The random noise generator was then used to put noiSe on the radar 
parameters, ATAE, EPSIE, RE, and RDOTE, and the apparent target location 
having coordinates XTE, YTE, and ZTE, was determined. Since a relatively 
small number of samples were used the mean, variance and independence 
of the noise for each of the parameters were checked using an autocorrelation 
routine. The results of analysis are shown in Tablel. 

The actual and noisy direction cosines for the three axé@S8' were then 
computed. These were, north direction cosine, DIRN, and DIRNE, east 
direction cosine, DIRE and DIREE, and down direction cosine, DIRD, and 
DIRDE. 

Knowing the probability of detection of the radar, RNDEV was used 
to determine a zero-one probability of detection for each frame time. This 
was done by cross plotting the absolute value of the deviatkhon required, 
versuS range, corresponding to the probability of detection curve, Fig. 2. 
This was then made piece wise linear and implemented by using IF statements. 
Thus for each actual range a maximum value of deviation for detection could 
be determined, if the absolute value of DEV obtained from RNDEV exceeded 


this value no return for that frame time was assumed. Fig. 6 illustrates the 


method used. 
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TABLE I 


AUTOCORRELATION ANALYSIS RESULTS OF PROGRAM 
THREE D 
(151 SAMPLES) 


eee Autocorrelation 


RN 1 

Target 

Acceleration 0) 10) LEY Sie 0.1787 
Noise 


RN 2 

Target 

Heading 

Acceleration | -0.001 1.0636 
Noise 


RN 3 
Target 
Pitch 
Acceleration -0.1032 0295015 
Noise 
ATA 
Antenna 
Azimuth 0.0785 eoG 12 
Noise 
EPSI 
Antenna 
Elevation 0.6308 1.0433 
Noise 
R 
Radar 
Range 0.1386 
Noise 


RDOT 

Radar 

Range 

Rate 0.1114 
Noise, ohn Se rate, 
a 
Radar @ © ey oe ee 
Probability 

of Detection -0.0564 
Noise 
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In the actual system this would be the status of the individual track 
in the TWS loop after being correlated and associated. 

A quality counter, NQ, was devised for the program which incremented 
in steps of two for each received return to a maximum of ten. For each 
missed return it degenerated by one. When the quality counter reached 
zero the track was dropped and no smoothing was performed until the next 
received return. 

The particulars of the smoothing and prediction section were as 
described earlier. 

The primary quality indicator for the system was the three dimensional 
rms smoothed difference from the actual position. Since high relative 
velocities between the missile and the target are considered very likely 
in the active phase of the missile flight it is necessary that the missile be 
positioned accurately at the beginning of this period so that its maneuvers 
may be kept to a minimum. In order to assist in achieving this a desired 
rms smoothing error limitation of approximately one-half mile inside 50 miles 
was imposed. On the output curves points were plotted indicating these 
limits. For range constant lines of points at plus or minus one-half mile 
were plotted. For azimuth and elevation points were plotted as a function 


of range uSing the relation: 
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Angular error Maximum = ARCTAN 


If all the parameters were smoothed to the desired value the three dimensional 
rms error, CIRSMO, could be as high as 0.875n.m. Thus on the curves of 
CIRSMO a line of points having this value were plotted. 

The program print out contained the run number, problem time, actual 
range, rms three dimensional actual error, CIRERR, rms three dimensional 
smoothed error, CIRSMO, a one or zero for detection or no detection, ND, 
and the status of the quality counter, NQ. When the quality counter was 
zero no contact was assumed and all nines were printed out. 

To observe the results of the smoothing method the sub-routine 
CALL DRAW was uSed to plot the interceptor and target positions in the 
north-east, N-E, plane, looking down from above, and the north-down, 
N-D, plane, looking in from the east axis. The subpoutine was also used 
to plot the difference between the actual and smoothed values of range 
and the direction cosines, which were plotted as angular differences, and 


the values of the three dimensional rms smoothed error. 











oe Results of the Analysis of the Proposed System 


Table II contains the important details of the problem profiles used 
in the analysis. The results of run 1A are shown in Table III and Figures 
7 to 15. The results of runs 2A and 3A are shown in Appendix D. 

Three basic profiles were considered. These were identified by the 
initial position and heading of the target with respect to the interceptor. 

For run 1 the target was positioned at 100 n.m. and was closing the 
interceptor. For run 2 the target was positioned at 10 n.m. and was 
moving away from the interceptor. In run 3 the target was positioned at 
approximately 60 n.m. ata 30 degree angle from the interceptor and was 
headed on a cross track with respect to it. The actual profiles flown by 
the target and the interceptor could be varied in a myriad of ways, however 
in reaching the conclusions that follow only those flight paths shown in 
Table II were used. 

The results obtained for the system, as proposed by Hughes, indicated 
that the smoothed tracks were still quite noisy. The direction cosine smoothing 
results obtained by analyzing the figures of these values compared to the 
desired limits indicated that the direction cosines were being smoothed inside 
the limits rather well. Some degradation was noted in the cross track cases 
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SMOOTHING RESULTS 
RUN 1A 


TsTIME, R NM=ACTUAL RANGE, CIRERR=SDIMENSIONAL RMS ACTUAL ERROR, 
CIRSMO=3 DIiMENSIUMAL RAS SMOGOTHEGSERRSR, NDO=b SETECTION, ND2e9 NO 
DETECTSON, NG=QUALITY“SCUNTER, WHEN “idd=9 TARGET CONSIDERED LOST AND 
ALL NINES ARE PRINTED, © °°" ‘°°? ae 


f 


ND N 


T R NM CIRERR CIRSMO Q 
. 100.5000 999.9999 999.9999 99 0 
2, 99,4871 999.9999 999.9999 99 0 
4, 98,4809 999.9999 999.9999 99 0 
6, 97,4805 2.8089 2.8089 ‘1 2 
8, 96,4819 5501 2,7924 oOo 14 
10, 95,4829 999,9999 999.9999 99 0 
12, 94,4898 999.9999 999.9999 99 0 
14, 93,4664 999.9999 999,9999 99 0 
16, 92,4470 999,9999 999,9999 99 0 
18, 91.4370 999.9999 999.9999 99 0 
20 90.9195 999.9999 999.9999 99 0 
22, 89,3964 999.9999 999,9999 99 0 
24, 88.3717 999.9999 999.9999 99 0 
26, 87,3597 999.9999 999.9999 99 0 
28, 86,3438 2.1180 2.1180 1 2 
30, 85 3354 1.6792 2.0969 0 1 
Gen 84,3289 999.9999 999.9999 99 0 
34, 83,3276 3.9493 3.9493 1 2 
36, 82 3209 2.5368 5.9162 0 4 
38, 81,3087 , 7825 2.4359 1 3 
40, 89,2970 32225 2.3294 0. 2 
42, 79,2848 9345 2.2242 0 fi 
44, 78,2667 999.9999 999.9999 99 0 
46, 77,2548 2.7971 2.7971 1 2 
48, 76,2567 2.4162 2,8147 0 1 
50. 75,2688 999.9999 999,9999 99 0 
52, 74,2774 999.9999 999,9999 99 0 
54, 73,2785 999.9999 999.9999 99 0 
56, 72,2765 7591 7391 { 2 
58, 71,2801 3.5995 1.1076 1 4 
60, 70,2918 2.9569 1.2629 0 5 
62, 69,2966 14,8983 1.4499 0 2 
64, 68,2999 2.9966 1.59692 1 4 
66, 67 3133 53755 7592 7 6 
68, 66,3347 3.4992 7694 0 5 
Fae 65 3566 2.4328 , 7920 0 - 
72, 64,3877 2.94979 1,6471 1 6 
74, 63,4255 4.4218 1.72635 0 5 
76, 62.4593 4.4005 1.8257 0 4 
78 61,4952 , 8885 1.9555 0 3 
80. 69,5331 ~8745 et 70 i 2 
82, 59 5727 3.2185 1.2012 0 4 
84, 58,6241 6.4126 22m. a Ou eS 
B6, 57,6929 1.9297 1.2447 0 2 
36 56.773 oe. C7 Vea we 0 ee al 
OT oe oo. 39> oe 2955 1.6465 > Ff 5 
JZ, "'S4,9099 1.0702 21,5794 - 9 2 
94, 53,9815 4.2629 1.7212 0 4 
96, 53.0494 999.9999 999.9999 99 0 
98, 52.1058 2.9109 2.9199 1 2 
100. Si 1811 2,3925 2,8871 0 i 
102, 99,2685 1.2764 2,1161 1 38 
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104. 
100. 
108, 
5 Be ae 
dines. 
114, 


118, 


20 
dee 
124, 
126, 
eee 
130, 
132, 
134, 
136, 
1358, 
140, 
142, 
144, 
146, 
148, 
130: 
152% 
154, 
156, 
124), 
169, 


162, 


178, 


182, 


188, 
190, 
192, 
194, 
196, 
198, 
200 
eye, 
204, 
206, 
2) ar, 


210. . 


21¢c, 
214, 
246, 
218, 
220s 
eae, 
224, 


49.5610 
48,4669 
47,5685 
46 6628 
45,7441 


14 4s, B idrd - 
4H, 87117 
42°, 9252" 


41,9624 
40,9986 
40,9442 
39,0998 
$8 ,1587 
37,2168 
36,2832 
557 5557 
34,4268 
33,4964 
32,5567 
31,6179 
30,6886 
29,7558 
28,8093 
27,8563 
26,8896 
25,9816 
25.1572 
24,4389 
23,8445 


23,3933 
23,0952 
22,9552 
22,9797 


(28,1515 


25,9222 


25,8224 


24.3980 
29,8556 
25.4395 
26,0370 
26,6251 
27,1899 
27,6986 
28,1555 
28,6017 
29,0550 
29,5092 
29,9700 
30,4558 
30,9019 
$1,3873852 
51,8487 
32,5278 


72.9686 


32,2933 


$u 7606 


54,2715 
34,7654 
35,2696 
35,7582 
$6 ,2585 


2.1280 
2.1478 
6511 
~7560 
6194 


oo 1, Cie 


1 27097’ 


a ek 


9878 
1.5656 
,9145 
8954 
.9466 
1.6227 
1.6066 
11,5899 
3722 
~ 4850 
6428 
1.8072 
1.19359 
1.4976 
1.1917 
8278 
,6589 
,8572 
2.5858 
»9185 
6228 


1.54904 
1.9082 
1,.799° 
2.4656 
1.8425 
1.8799 
2 , O12 7 

2g? 
1.1249 
2,23579 
1.1206 
1.1881 
1.2085 
1.8540 
2.0981 
1.8206 

,8297 

~9816 

. 3614 

4553 
1.9997 
1,2539 
1.2489 
2,3558 


. 2.6509 - 
“1. 0851 


1.1820 
1001 
1.1220 
, 3621 
7449 
1.0007 
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246, 


222, 
254 
256 
2D Be 
269, 
262, 
264, 
266, 
268, 
270. 
272, 
274, 
276, 
278, 
eon 
282, 
284, 
286, 
288, 
290, 
292, 
294, 
296, 
298, 
S00, 


¢ 


36.7614 
37,2662 
37,7737 


gma, 28220, 
| 38,7943, 
397369 > * 


$9 ,8210 
40,5578 


40,8568 
‘41,3762 


‘41,8975 


42,4293 
Maen 
43,4794 
43,9969 
44 5255 
penne 
45 5847 
46,1162 
46 6486 
47-1547 
47,7164 
98,2528 


~48 7889 


49,3267 
49,8658 
50,4093 
90,9496 
91,4871 
oe ege 
D2 os ce 
2S i130 
53,6600 
54,2952 
54,7504 
55,2968 
25 8449 
56,5915 


TABLE III 

2.1929 ~6306 
4295 .9070 
eed ,8871 
2.2299 ,e91k 
MA PAOT aay OG 
SV272°" > FTC 
,9596 ,208S89 
, 8785 ,7262 
2,/7186 1.19638 
6,0221 1.2214 
1.1927 1.2999 
$,5265 1.2827 
,9221 1.0891 
22,2875 25461 
1.9963 ,2977 
32446 tS 272 
4.0991 ee 10 
4.3779 14.2842 
S$ ,1691 1.8915 
2.8110 2.2411 
5.2857 .9999 
5.9789 2,6170 
1,S971 179559" 
3.9975 , 7339 
35,3489 .8522 
4459 .4249 
1.0696 6642 
$,8562 2.0458 
4.9476 1.1952 
,8078 ~5816 
1.3545 ~63596 
1.9762 9642 
4.8418 1.0683 
1,8616 1.1776 
2,180 41,2917 
1,1498 ~6087 
2.4983 ~6028 
41,6145 6383 
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X-SCALE = 5.Q2E+814 UNITS/INCH, 
¥-OCALE = 4J.QQ6-22 UNITS7 INCH 


GENT2 EAST SMOOTHING RESLILTS. © RUN 1A 
SANRIDR NAIA 
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FIGURE 14 





XSCALE = S,QQE+21 UNITS INCH, 
“-ICALE = 1,QQE-22 UNITS INCH. 


GENTZ BORN SMOSTHING- RESULTS RUN 14 
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A-OCALE = 5.QQ86+81 UNITS INCH. 
4M-SCALE = 6.QQE-84 UNITS/ INCH. 


CENTZ. 30, RMS, SMOOTHING RESULTS RUN 1A 
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where the rates of change of the direction cosines were high. As was 
expected the range smoothing results were the prime cause of the smoothed 
error. 

The system was able to smooth the tracks to within one and a half 
nautica! miles inside 50 miles rather consistently, but it tended to respond 
rather strongly to the noise. It was noted, however, that this had the 
effect of causing no degradation in performance when the target began 
maneuvers. Also all the profiles yielded very similar smoothing results 
indicating that the smoothing constant values uSed were near optimum for 
the fixed constant, non-adaptive, system. 

In an attempt to more closely define the optimum values of the 
smoothing constants a program which varied the constants for each smoothed 
parameter from 0.05 to 1.00 in steps of 0.05 was written. This program 
is named ALPHA CHECK and listed in Appendix C. The run profiles used 
were Similar to those listed earlier, but, the run time was shortened to 
100 seconds so that only the time that the target maneuvered was considered. 
The average mean square error for the complete run was computed for each 
smoothed quantity over the range of the values of its smoothing constant. 

An autocorrelation analysis was also performed to evaluate this noise 


distribution.  Thevorosults ara shoyvm in>,Table IV... Eor these runs the 
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AUTOCORRELATION ANALYSIS RESULTS OF PROGRAM 
ALPHA CHECK 
(51 SAMPLES) 


Variance Largest Value of 
| Autocorrelation | 
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the probability of detection was one. This did not effect the results since 
no smoothing would have been performed when a contact was missed. The 
output would have been affected in an indeterminable way, however, by 
these predicted values. 

The results of these runs are shown in Table V. It was apparent 
from these results that there was considerable interplay between all the 
variables. Some rationale could be seen in this interaction, but it appeared 
nonlinear, since any change in the target's profile varied all the points of 
the optimum values for the smoothing constants. It was noted further that 
the knuckle in the curves at the optimum point was quite sharp indicating 
that system performance could be seriously affected when improper values 
of the smoothing constants were used. 

Due to the variations which were found in the optimum values of the 
smoothing constants for the chosen profiles it was not believed to be 
necessary to vary the parameters more widely since it was apparent from 
these runs that the system was very sensitive to changes in the profiles 
or the values of the system noise. 

The question must then be asked as to whether the reQuirements 


established for the analysis were too strict. In this regard specifically 
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the standard deviation of range noise and the 0.5 mile value chosen for 
the desired smoothing error become suspect since they were somewhat 
arbitrarily chosen. In rebuttal itis strongly felt, in the light of the high 
overall performance required of the entire system, its cost, and purpose, 
that the use of this relatively uncontrollable method of smoothing could be 
disastrous. Since the output of the smoothing section has considerable 
control over the performance of the system, Specifically the accurate 
positioning of the missile, it is believed that a more adaptive system is 
desirable. It is further felt that a programmed logic system which would 
provide for varying the values of the constants would not be warranted 
since more sophisticated, and proven adaptive systems are-available. 
Possible choices for such a method could be either the Least Squares, 
or Kalman filter methods. 

In short it must be concluded that the “yv, @ ; % Tracker 
system does not have the response or versatility required for a high 


performance air to air intercept problem. 


>’ >. 


wwe @ 
ed 
w 


% 

} & 

> a) 9 
9 





43 of 100 


aw 


ccc 


¢ 


eecee 


qeea 4 





eeee 4 


¢ 8 + a ® & 6 8 aes %@ ee @¢ } e¢ 
, Oo ¢ 5 a i 
7 2 fi y 
e°¢ ee 


F . : e? e? , ee) 
eee ° 


BIBLIOGRAPHY 


Hughes Aircraft Co., AWG-9 System Evaluation and Analysis 
Report, June 1964. 


Skolnik, M.I1., Introduction to Radar Systems. McGraw-Hill 
Book Co., 1962. 


Benedict, T. R and Bordner, G. W. Synthesis of an Optimal Set 
of Radar Track-While-Scan Smoothing Equations. Institute of 
Radio Engineers Transactions on Automatic Control. 

Vol. AC-%, 27-32, July, 1962. 


Simpson, H. R. Performance Measures and Optimization Condition 
for a Third-Order Sampled-Data Tracker. Institute of Radio Engineers 
Transactions on Automatic Control. Vol. AC- , 182-183, April, 1963. 


Skalansky, J. Optimizing the Dynamic Parameters of a Track-While- 
Scan System. RCA Review. Vol.XVIJI, 163-183,June, 1957. 


; jm , 
®a)>e0e ) } ® » 6 e * } ® ° 6 a 8 e GaP POL N TTA 
e ) 
) 


e 
9 ° ‘ ) @ ) 
6 e ) ® 6 ) 6 
® 


5 EI so com 


® 
® ® m ¢ (a 
4 





4 


4 


eee 4 


« 4 


4 


4 


« 


eC i 


a ¢ 


teed 4 


io 





@eo6 


NOTE: 


ADEV: 


ADSMO: 


AESMO: 


ANSMO: 


AE: 


AN: 


AERRM: 


AERRP: 


APPENDIX A 
TABLE OF SYMBOLS USED IN COMPUTER PROGRAMS 


If a symbol is used only in the ALPHA CHECK program it is so 


noted in parenthesis. 


Absolute value of deviation obtained from RNDEV. 


This is uSed in the probability of detection computations. 


Down, east, and north direction cosine differences 


from actual, and converted to radians. 


(ALPHA CHECK) Values of alphas for down, east, and 
north direction cosines, and for range, uSed for determining 


the optimum values. 


Negative and positive desired limits on ADSMO, AESMO, 


and ANSMO, a function of actual range. 
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ALPHAD: 
Pee. 
ALPHAN: 


ALP HAR: 


ANOISE: 


AT: 


ATDOT: 


ATA: 
-ATAE: . 
ATAD: 


ATADE: 


BETAD: 
BETAE:? 


BETAN: 


Bae 


Bape: 


Position smoothing constants for down, east and north 


direction cosines, and for range. 


(ALPHA CHECK) Relative amount of radar parameter noise 


used in problem. 


Target heading and turn rate used in computation of target's 


position. 


Antennas actual and noisy azimuth angle meaSured positive 


clockwise, in radians and degrees. 


Rate smoothing constants for down, east, and north 


direction cosines. 


TargetS pitch angle and pitch angle rate used in computation 


of targets position. 
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CIRERR: 


CIRSMO: 


DEVE 


DT: 


Otis” 


DIRNE: 


DIRNS: 


DIRNP: 


ey ih tess 


PIREE: 


DIRES: 


DEREP: 


IRD: 


ITRDE: 


DIRDS: 


PIRIOL: 


Three dimensional RMS actual and smoothed error difference 


from actual target position. 
Random deviation obtained from random number generator. 
Target position computation step size. 


North direction cosine actual, noisy, smoothed and 


predicted. 


East direction cosine actual, noisy, smoothed and 


predicted. 


Down direction cosine actual, noiSy, smoothed and 


predicted. 
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DIRND: 


DIRNDD: 


DIRED: 


DIREDD: 


DIRDD: 


DE.DD: 


DDSMO: 


DESMO: 


DNSMO: 


DELTAT: 


DIVE: 


EPSID: 


BE eQHDE: 


EPol: 


BPOEE: 
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North direction cosine, rate and acceleration. 
East direction cosine, rate and acceleration. 
Down direction cosine, rate and acceleration. 


(ALPHA CHECK) Down east and north direction cosine 


smoothing differences from actual. 


Target decceleration time to final velocity at 1.75 g. 
Target maximum dive angle, radians, corresponds to 60 deg. 


Antenna actual and noisy elevation angle measured positive 


downward, in radians and degrees. 
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ERAD: 


ERAE: 


ERAN: 


ERAR: 


FPOINTS: 


GAMMAD: 


GAM MAE: 


GAMMAN: 


rr: 


K: 


(ALPHA CHECK) Average mean square smoothed error for 
down, east, and north direction cosines, and for range. 
They are the quality indicators used for obtaining the 
optimum values of the smoothing constants. 

A counter which fixes the target's pull out velocity. 
(ALPHA CHECK) Number of radar returns received. 


4 


Acceleration smoothing constants for the down, east, 


and north direction cosines. 


Target maneuver load factors for push over, pull out, 


and turn. Note: GDN is a negative number. 


Hollerith characters which identify the output plots. 


Number of frame times since problem start. 


Number of desired limits points plotted from problem start. 


* Read into the program 
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ND: 


NERMAX: 


N GRAF 1:* 
NGRAF2:* 


N GRAF 3:* 


NOISE:* 


N PRINT: 


NQ: 


NQC: 


) ’ ’ 

s } ’ 5 ’ ? ne) 
’ ) 
) 


(ALPHA CHECK) Number of values of smoothing constants 


used in finding optimum value. 


A zero or one depending on whether detection was missed 


Or not. 


Number of points plotted indicating the desired limits on 


the smoothed errors. 


Setting NGRAFS to zero supresses various portions of the 


graph output. 


Number of times RNDEV is stepped prior to problem start to 


allow a different sequence of random numbers to be called. 
Number which separates pages of print out. 
Quality counter. 


Counter which reinitializes the initial conditions after the 


quality counter goes to zero. 


* Read into the program 
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NUMPTS: 


NUNIF: 


Onl 


ED: 


A: 


POVERD: 


POUTD: 


PSIO: 


PSIOD:* 


PSIT: 


Bol. * 


PSITO: 


@: 


* Read into the program ,, 
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Number of curve points plotted. 
Counter in random number generator. 
Intercept:r Mach number. 


Target detection probability for each frame time, a function 


of the actual range. 

S214 159265. 

Target pushover rate, rad/sec. 
Target pull out rate, rad/sec. 


Interceptor heading in radians and degrees. 


Target heading. 


Target initial heading, deg. and rad. Note a zero not an 


alphebetic o on target’s initial positions. 


Number of frame,times since last radar return of range rate. 
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Re Actual and noisy range to target, ft andn.m. 


Es: 

RN M: 

RN ME: 

RERRM: Negative and positive desired limits on range, a 
eeESINE constant .5n.m. 

RDOT: Actual and noisy range rate of target with respect to 
RDOTE: interceptor, positive opening, ft/sec. 

RDOTM: Range rate actual, noisy, predicted, and smoothed, 
RDOTME: nome, sece 

RDOT MP: 

RDOTMS: 

RDDOTE: Target noisy acceleration, atm. QE" ; 

RN1: Random accelerations of target for velocity, heading, 
RN2: ana piten. 

RS: 
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RN MP: Predicted and smoothed range to target, n.m. 
RN MS: 

RSMO: Smoothed range error from actual, n.m. 

ae Time since problem start. 

TA: Total target iteration time since problem start. 
hCG: Iteration time which determines time for plotting 


desired limit points. 


2 Time to target dive. 
Tvl: * Target Mach number. 
iS: Check time against iteration time to enter main program 


from target iteration. 


aes Time to target turn. 

vo: Intercepter velocity, ft/sec. 

va: Target velocity and acceleration, ft/sec. and ft /sec? ; 

Viper: 

Vie : Final target velocity, a function of the target's final 
altitude. 

* Read into the program 
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XONM: 


XTERR: 


xLO:* 


XTE: 


XTONM: 


XTNME: 


XTNMS: 


XTSMO: 


Nel: 


Ze 


Points which mark initial target location on the smoothed 


plot. 


Interceptor position along the north axis, ft. and n.m. 


Difference between target's noisy and smoothed position, 


nem. 


Target position along north axis, initial, actual, and noisy 


in feet, and initial, noisy, and smoothed position, n.m. 


Target's difference between actual and smoothed x position. 


Similar to Xl and X2. 


* Read into the program 
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YTERR: 


WO :* 


ole: 


YTONM: 


YTN ME: 


YTNMS: 


YTSMO: 


ZAMS 


L2: 


ALLA: 


ZO °* 


ZONM: 


ZTERR: 


C 


Similar to XTERR. 


Target's position along the east axis, similar to X positions. 


Target's difference between actual Y position and smoothed 


Y position. 


Similar to X1 and X2. 


Target's deflection angle from the horizontal, measured 


positive downward from the Intercepter, rad. 


Interceptor's altitude, ft. andn.m. 


Similar to XTERR. 


* Read into the program 
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ZTO:* 
aS 
AE: 


ZINM: 


ZINME: 


ZINMS: 


ZISMO: 


Zee 


Ze: 


Target's altitude similar to X positions. 


Target's difference between actual Z position and smoothed 


Z position. 


Final target altitude, ft. 


Target pull out altitude, ft. 


* Read into the program 
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APPENDIX B 
COMPUTER PROGRAM DETAILS 


A. Purpose: The program is designed to generate three dimensional 
target-interceptor geometries for use in analyzing optimum smoothing 


schemes. 


B. Limitations: The interceptor starts at the coordinate axis origin and 
may have any northerly heading. The target may be placed anywhere in 
the northern semi-circle with any initial heading. The positive axis 
directions are: X, north, Y, east, and Z, down, measured from the 
interceptcr's altitude. The problem will stop anytime the target flies 


out of a 65 deg. cone from the interceptor's heading. 


C. Program Capabilities: The program will compute the actual range, 

range-rate, antenna azimuth and elevation for a radar frame time of 

two seconds. It will put noise on the above having the standard 

deviations of 2.5 n.m. for range, 8.72kts for range-rate, 0.745 deg. 

for azimuth, and 0.86 deg. for elevation. It also puts a half g acceleration 

noise on the target. It will compute the noipy target location, based on the 
a ee 
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noise. It will further transfer the actual and noisy radar parameters 

into a moving north, east, down coordinate system which is fixed to 
move with the interceptor. It will determine whether detection was made 
or not each frame time, based on a detection probability curve, and 
perform a quality counter operation which prevents non-detected targets 


from being smoothed after they are lost for a reasonable number of frames. 


D. Input data requirements: Four input cards are required. 

1. The first card reads in 16 Hollerith characters which identify the 
program on the print out and graphs. The format is 2A8. 

2. The second card reads in the interceptor's altitude in feet, the 
target's X and Y coordinates in nautical miles, and altitude in feet, 
and interceptor and target heading in degrees. The format is 7F6.0. 

3. The third card reads in the target and the interceptor’s Mach number, 
the target's maneuver g load capabilities for turn, pull out, and dive, 
(negative number), and the problem time to target turn and/or dive. 

The format is 7F6.2. 

4. The fourth card reads in a fixed point number that allows a different 
Starting point for the subroutine RNDEV, and three fixed point ones or 


zeroes which control the graph output. The first, if zero, suppresses 
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the profile curves, the second, if zero, suppresses the smoothed 
difference plots,,and the third, if:zero, Suppresses the three dimensional 


yms smoothed difference plot. The format is 416. 


E. Requirements for reprogramming the smoothing section: The entering 
arguments of the smoothing and prediction section are RNME, DIRNE, 
DIREE, DIRDE, and NOC. The suggested statement numbers to uSe are 


the 600 series or above 2000. 


F. Quality indicators: The primary quality indicator is the three 
dimension yms smoothed error from the actual position. The three 
dimensional rms noise error is computed for comparison. The desired 
smoothing error on the individual quantities is 0.5 n.m. This could 
give a three dimensional rms smoothed error of as much as 0.875 n.m. 


The output curves have points plotted indicating these limits. 


G. Output: The print out contains the run identification read in, the 
problem time, the three dimensional Y¥Yms actual and smoothed error, 
whether detection was made or not, ND equals one or zero, and the 

- Status of the quality counter, NQ. There are nine output curves. They 


plot the following: 
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1. The interceptor and target's actual and noisy positions in the 
north-east plane. 

2. The interceptor position and the target’s smoothed position in the 
north-east plane. 

3. The interceptor and target’s actual and noisy positions in the north- 
down plane. 

4. The interceptor'’s position and the smoothed target's position in the 
north-down plane. 

5. The range smoothed differences from actual and the desired limit 
points, versus time. 

6. The north direction cosine smoothed differences from actual and the 
desired limits points plotted as a function of actual range, versus time. 
7. The east direction cosine smoothed differences similar to the north 
direction cosine curve. 

8. The down direction cosine smoothed differences similar to the north 
direction cosine curve. 

9. The three dimensional yms smoothed error from actual and the 


desired limits points, versus time. 


H. Runtime: The average run time is between 2 minutes and 30 seconds 
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PROGRAM THREE D 


THIS PROGRAM COMPUTES THREE DIMENSIONAL INTERCEPTER-TARGET 
GEVMETRIES FOR THE PURPOSE’ OF) OP TIMTAING: VAROUS PARAMETERS, 

THE PROGRAM WILL BRP ANY” TIME THAT ‘THE TARGET FLIES OUTSIDE OF A 
65 DEGREE GONE’ FROM ‘THE INTERCEPTERS ' HEADING. 


FOUR INPUT CARDS ARE REQUIRED 

bie THE FIRST CARD READS IN 16 HOLLERITH GHARACTERS WHICH IDENT= 
IFY THE PROGRAM ON THE PRINT OUT AND GRAPHS. THE FORMAT IS 2A8, 
ee THE SECOND CARD READS IN INTERCEPTERS ALTITUDE IN FEET, 
TARGETS X AND Y COORDINATES IN NAUTICAL MILES AND Z COORDINATE 
IN FEET, AND INTERCEPTER AND TARGET HEADING IN DEGREES. THE 
FORMAT IS 7F6.0. 

3. THE THIRD CARD READS IN TARGET INTERCEPTER MACH NUMBER, 
TARGET MANEUVER G LOAD CAPABILITIES FOR TURN PULL OUT AND DIVE 
(NEGATIVE NUMBER), AND THE PROBLEM TIME TO TARGET TURN AND DIVE, 
THE FORMAT IS 7F 6.2. 

4. THE FOURTH CARD READS IN A FIXED POINT NUMBER THAT ALLOWS 

A DIFFERENT STARTING POINT FOR THE SUBROUTINE RNDEV, AND THREE 
FIXED POINT ONES OR ZEROS WHICH CONTROL THE GRAPH OUTPUT. THE 
FIRST IF ZERO SUPPRESSES THE PROFILE CURVES, THE SECOND IF ZERO 
SUPPRESSES THE SMOOTHED DIFFERENCE PLOTS, AND THE THIRD IF ZERO 
SUPPRESSES THE 3D RMS SMOOTHED ERROR PLOT. THE FORMAT IS 416, 


DIMENSION AND READ STATEMENTS 


DIMENSION 17(451),XONM(151)>YONM(151),>,Z0NM(2151)5XTNM(151),YTNM(151) 
1,Z7NM(6151)+RNM(151),ATAD(151),RNME(151),RDOTME(151),ATADE(151), 
2RNMS(454),RD0TM(454),XTNME(151)5YTNME(154),ZTNME(4151)5XTNMS(14514), 
SYTNMS (151),ZTNMS(1591), 17(12),DIRN(191),DIRE(151),DIRD(151), 
4DIRNE(151),DIREE(191),DIRDE(151),RDDOTE(191),RNMP(1924),RDOTMP(151) , 
5,DIRNS(154),DIRES(151),DIRDS(191),DIRND(151),DIRED(151),DIRDD(151) 
6,DIRNDD(154),DIREDD(151),DIRDDD(151),D0I1RNP(151),DI]REP(191),DIRDP 
7(151),EPSID(451),EPSIDE(151),RSM0(151),ANSM0(1514),AESM0(454),ADSMO 
8(454),ZETA(151).NQ(154),XTSMO0(191),YTSMO(191),ZTSMO(191),CIRERR(15 
94),AERRP(39),AERRM(59),RERRP(39),RERRM(5Q),TC(59),RERR(191) 

DIMENSION X1(2),Y1(62),21(2),X2(2),Y2(2),Z2(2),CERRP (30) 
1»>XTERR(454),YTERR(154),ZTERR(151),CIRSM0(151),ND(1514) 

DO 2 J=1,6 

READ 100¢(IT(I),J]=5,6) 

100 FORMAT(2A8) 

READ 11,>Z0,XTOQONM,YTONM,ZT9,ZTF,PSIOD,PSITD,™TM,0OM,GT,GDP,GDN, 

17T,TD,NOISE,NGRAFA,NGRAF 2, NGRAFS ; 
11 FORMAT(7F6.0+/4,7F 6120/0416) 


INITIAL CONDITIONS 


500 T(1)=0, 

TC(1)=0, 

TARO + Gy ee ae poe ee 
al) i ee 7: oe : 
NERMAXESQ 78> 28 ee weer ree bk eee 5 Chae 
P1=3,14159265 . 

RAD=180./PI 

DIVE=-60./RAD 

DT=.09 So 
NUMPTS2454 


NPRINT=68 _ ed 





ccet ¢ 


c 


t 





Sd 
eo 
eo 
= 


63 


69 
57 
73 


74 
75 


78 


719 


58 


122 


102 
109 


22 


Rie a tae y(t 
XTO=XTONM#6076, 4 - 
YTO=YTONM¥6076,1 

XTsXTQ 

YT=YTO 

LWEZ10 

PSITQ=PSITD/RAD 
PSIO=PSIOD/RAD 

PSIT=PSITQ 

AT=PSITO 

ATDOT=0, 

ATDDOT=0. 

BTl=0. 

OG Ue 

BTDDOT=0. 

RNL=Q- 

RN2=0, 

RNS=Q. 

NUNIF 212297935125 


TARGET AND INTERCEPTER VELOCITY COMPUTATIONS 


Varin) pave 
IF (ZT0737000. )65,69, 69 
VT=TM*(1116.-.00598"7TQ) 
GO TO 57 
VT=TM*967,6 
1F(Z20-3S7000.)75,74,74 
VO=OM*(1116.-,.00398"Z0) 
G0 10 75 
VO=0M*967.6 
IF (ZTF-10000.)78,79,79 
VTF=1340.+0-088ZTF 
DELTAT=(VT+VTF)/96,4 
GO To 58 
Vie=V I ; 
DELTAT=0Q. | 
CONTINUE 


THIS ALLOWS A VARIATION IN THE: RANDOM NUMBER SEQUENCE 
[F(NOISE)199,109,122 “ 

DO 492 L=4,NOISE 

CALL RNDEV(NUNIF, DEV) 

CONTINUE 


BEGIN MAIN PROGRAM DO. LOOP 
DO. 2g Lesea Sar ee pe at 
CHECK’ I'VERATION TINE AGATNST RADAR FRAME TIME 
[F (ABSF(TS-TA)-10,E~-9)29,29,24 
CHECK MANEUVER TIMES AND/OR READ IN MANEUVER CONDITIONS 


62 of 100 


anee ‘ 


e888 & 


ee 


ee ae 





eeee cae 


aeae 





24 
39 


oo 
29 
40 


a4 
76 


77 


se 
46 


45 
54 


35 
2038 


37 
200 


204 
43 


38 
205 
206 


202 


24 


22 


59 


TA=TA*DT 

[FC TA-TT)39, 28,286 

ATDOT=0. 

GO TO 29 

ATDOT=GT¥32.47/VT 

TFTA? 7D 749 240,44 oe ame ee ee ee ee eee 
PO ee : re on ee va : or . 
G0 TD Roe prrvy 0 » YOO yO yr? 9 ve: ae cae 


IF CZT-ZTF)76,77,77 
BTDOT=(GDP-1.)*32.,47/VT 
GO To 42 
BIDOT=(GDN-1.)%32,147/VT 
VTDOT=-56,4 

IF (CCPSITQ+*#PI)-AT)496,496,45 
oon a0 

AT=PSITo+PI 

IF (BT-DIVE)34,34,35 
a0 1 = 0% 

BT=DIVE 

IF (F) 201,291,208 
ZTPO=(1,-COSF(BT)) *VTee2/((GDP~1,.)*32617)+ZTF 
IF (ZTPO-ZT)38,37,387 

IF (F 204,201,200 

FoF-4. 

VTPO=VT 
BTDOT=(GDP-1.)*32,47/VTPO 
IF (87T)38,438,438 

BIDOT=9,. 

BT=0. 


IF (ZT-ZTF)205,205,206 

ZT=ZTF 

IF (TA-(TD*DELTAT) )91,51,292 7 
VTDOT=0. 

VTSVIF 


ITERATION OF TARGETS POSITION 


VT=VT+VTDOT#*DT+RNi «DT 
ATZ=AT+ATDOT*¥DT+RN2eDT ev 2/2, 
PSIT=AT 
BT=BT+BTDOT*DT+#RNSeDT «x2/2, 


COMPUTE TARGET POSITION 


XT=XT#VT*COSF (BT) ®COSF (AT) *DT 
YT=YT+#+VT«*COSF (BT) *SINF (AT) #DT 
ZT=ZT+VT*SINE (BT) wDT 


RETURN TQ START TO CHECK AGAINST RADAR LOOK TIME 


GO TO 22 
TS=TS+t2, 


COMPUTE. .INTERCEPTERS, POSITION, AND #&GTUA! RADAR, PARAMETERS 


» 700 4 apo 9 9” > o v ) ’ o> 94 » 
¢ ’ ’ ) ) ) » > ’ ) 


XOSVOR FCI #COSFCPSICH: ore re) Oo, ee 
YO=VO*T(1)*SINF(PSIO) 
R=SQRTFCCXT"X0) wwe ( YT HY) wee 2e( ZT 920) #82) 
EPSI=ASINF((Z0-ZT)/R) 

IF (PSIO=P1)69,59,59 

PSIO=PSJ]0-2.«P] 
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60 


32 


$3 
26 
49 
48 
29 
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14 





IF CYT-Y0)51532,52 

ATA=-ACOSF (CXT-XO)/(R*COSF(CEPSI)))-PSIO 
GO 10 33 

ATA=ACOSF ((XT=X0)/(ReCOSF (CEPSI)))-PSI0 


YG Pe" TF: TARGE cs) STILE: TY RADAR: BEAM, IF NOT STOP 


y ,’_v 3 >? 


[FCRTACL AS4Sj26,BOISg cet ta ta 
IF CATA*2.13495)50,49,49 

IF CEPSI=-1,1345)48,48,59 
IFCEPSI+1.1545)50,27,27 

NUMPTS=([-4) 

NERMAX=NUMPTS/19 

NPLOT=NUMPTS-NCTC 

IF (NUMPTS~68)19,14,14 

NPRINT=NUMPTS 

GO TO 1 


27 RDOTH-VT*COSF (PSIT@ACOSF ((XT=X0)/(R*COSF(EPSI)))+PI)*COSFC(EPSI) 


70 


1-VOxCOSF CATA) «COSF (EPS!) 


CONVERT ACTUAL RADAR PARAMETERS TO KTS, NM/SEC, DEG 


RNM(I)=R/6076.14 
ATAD(I) SATA*RAD 
EPSID(I)=EPSI*RAD 
RDOTM(])=RDOT/s6076,1 
XONM(1)=X0/6976,1 
YONM(1)=Y0/6976,1 
ZONM(1)=Z0/6076,1 
XTNM(1)=XT/6076,1 
YTNM(1)2=YT/6076.1 
ZTNM(1)2ZT/6976,1 


OBTAIN TARGET MANEUVER NOISE 


CALL RNDEV(NUNIF,DEV) 
RN1=16.*DEV 

CALL RNDEV(NUNJF, DEV) 
RN2=15./VT*DEV 

CALL RNDEV(NUNJF,DEV) 
RNS=16,/VT*DEV 


OBTAIN NOJSY RADAR PARAMETERS AND CONVERT TO CONVENIENT UNITS 


CALL RNDEV(NUNIF, DEV) 

ATAE=ATA+,Q13*DEV 

CALL RNDEV(NUNIF, DEV) i 
RE=R+15990-*DEV , 
CALL RNDEVC(NUNIF, DEV) 

EPSIE=EPSI+.919«DEV 

CALL RNDEV(NUNIF, DEV) 

RDOTE=RDOT+14.7174 «DEV 

ATADE(I)=ATAE*RAD 

EPSIDE(1)=EPSIE*RAD 

RNME(1)=RE76976,4 ops Vee) ae 
RDOTME(I:) = RDITE £69 7644 : a0 0 0 8 » 0 


COMPUTE NOISY TARGET SPACE STABILIZED POSITION 


IF(PSIO+1,)79,70072 


XTE=XO+RE*COSF CATAE+PSI0) #COSF CEPSIE) 
YTE=YO*REwSINE (ATAE+PSI0) «COSF (CEPSIE) 
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66 
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68 
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53 


250 
54 


254 


Somos) 2 

XTE=XO+RE*COSF (ATAE=PS10) #COSF (EPSIE) 
YTE=YO+RE*SINF (ATAE-PSIO) #COSF (EPSIE) 
ZTE=ZO-RE*SINF (EPSIE) 

XTNME( 1) =XTE/6076.1 

YONMER 1 Voice (0p 7 Oude) cen 
ZINMEXP=ZTEA6O7G, 99>, > tds fe 2 


y 
é y } y 
? $9994 y d eovy ,vV°9 Yeo e ’ a 





via) 


COMPUTE ACTUAL DIRECTION COSINES FOR ‘INTERCEPTER WITH ROLL 
AND PITCH ZERO 


DIRN(I)=COSF (CATA) ®COSF (PSIO)=SINF CATA) & SINE (CPSIO) 
DIRECIL)=COSF CATA) ®SINFE(PSIO)*SINE CATA) #COSF (CPSIO) 
DIRD(I1)=-SINFCEPST) 


COMPUTE NOISY DIRECTION COSINES 


DIRNECI)=COSF (ATAE) *COSF (PSIO)-SINFCATAE) #®SINF (PSIO) 


DIREE(CI)=COSFCATAE) *SINF (PSIQ)+SINF (CATAE) *COSF (PSJO) 
DIRDECI)="-SINFCEPSIJE) 


OBTAIN DETECTION PROBABILITY AND DETERMINE Jb TARGET WAS 
DETECTED 


CALL RNDEV (NUNIF,DEV) 
IF(RNM(1)-20.)21,25, 23 \ 
PD=2,36-.043*RNM(]) 

GOTO 52 

IF (RNM(1)*99,)36,47,47 
PD=1,95~.0239*#RNM(1) 

Go To. 52 

IF (RNM(1)=89.)65,66, 66 
PD=1.55~.016*RNM(J) 

GO TO 52 

IF (RNM(1)-110,)62, 668,68 
PD=.88-.0077*RNM(1) 

GO TO 52 
PD=.25-.0019*RNM(I) 
CONTINUE 

ADEV=ABSF (DEV) 

IF (PD~ADEV)53,54,54 


QUALITY COUNTER COMPUTATIONS 


ND(1)=0 

NQ(1)=NQ(I71)=2 

IF (NQ(1))259,151,152 . 
NQ(1)=0 

GO To 153 

ND(1)=1 

NQ(1)=NQ(J-1)+2 

IF (10-NQ¢1) 254,153,153 

NQ(1)=10 

GO TO 153 


»>)>) 7.00 


ae SMOOTHING PND PRENICTION SECTION: 
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Own nent > 3” Satses 5 a7” % i , ) 


NO DETECTION CONDITIONS 


IF QUALITY COUNTER IS ZERO PRINT THE FOLLOWING AND START NEXT 
FRAME TIME 





ww 


eed 


{ ce (cere 


€ 


€ 


eww 





‘ 
‘ 
€¢ 
‘ 
‘ 


ew@wvde 


oto 


wea, « 


weve@ed 





14 


152 
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154 


64 


CIRERR(1)=999.9999 

CIRSMO(1)=999,9999 

XTNMSCT)=XTNM(J) 

YTNMSCL)=YTNMCT) 

ZTNMS(T)=Z1NM(T) 

NDCIIE99 i) 0° 2) 2 a ee 
NQC=1: — ald ee ewe ot NS 
Go To 20° 20 oe ee) 


TF QUALITY COUNTER IS NON ZERO ERROR VALUES EQUAL PREDICTED 
ONES, COMPUTE ONLY NEXT PREDICTED VALUES 


f 
RNME(1)= RNMP (1) 
RDOTME(1)=RDOTMP(]) 
DIRNE(I)=DIRNP(J) 
DIREECI)=DIREP(1) 
DIRDE(J)=DIRDP(]J) 
GO TO 64 


DETECTION CONDITIONS 
IF (NQC)64,64,154 


INITIAL CONTACT CONDITIONS AT START OR AFTER CONTACT HAS 
BEEN LOST AND REGAINED 

NQC=NQC-1 

RNMP(I)=RNME(T) 

RDOTMP(I)=RDOTME(CT) 

DIRNPCJ)=DIRNE(CT) 

DIREPCI)=DIREEC]) 

DIRDPCI)=DIRDEC]) 

RDDOTE(I-1)=0. : 

DIRND(I-1)=0. 

DIRNDD(I-1)=0. 

DIRED(I-1)=0. 

DIREDD(I-1)=0. 

DIRDD(I-1)=0. 

DIRDDD(I-1)=0. 

CONTINUE. 


SMOOTHING CONSTANT COMPUTATIONS 


ALPHAR=,49 

ALPHAN=.43 

BETAN=ALPHAN#*®*2/(2.-ALPHAN) 

GAMMAN=2. *BETAN/ALPHAN=(ALPHAN+BETAN ) 
ALPHAE=ALPHAN 


ALPHAD=.3$5 
BETAD=ALPHAD**2/(2,-ALPHAD) 


GAMMAD=2. *BETAD/ALPHAD-(ALPHAD+BETAD ) 
SMOOTHING AND PREDICTION COMPUTATIONS 


nena 2 290 Oo nO >) i 
RNMS(1)= RNMPC LeCAL PHAR S GHNMBLT3 RNMPCLDD. oo 
RDDOTEC(I)=RDDOTE( I -1)41./(2. 0) *(RDOTME( J )&= “RDOTMP(J)) 
RNMP(1+1)=RNMS( 1) +2, *#RDOTME(])+2,*RDDOTE(]) 
RDOTMP(1+4)=RDOTME(1)+2,*RDDOTE(]) 


NORTH DIRECTION COSINE 
OKIE 
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210 
212 


ead 
220 
222 


221 
250 
2352 


254 


DIRNSC(CT)=DIRNP(])+ALPHAN®* (DIRNECT)-=DIRNP(I)) =, 


DIRNDCT)SDIRND(1-1)+2.*DIRNDD(1-1)+BETAN/2.,*(DIRNE(I)-DIRNP(I)) 


DIRNDDCI)=DIRNDD(1-+1)+2.*GAMMAN/4.,«*#(DIRNECI)-DIRNP(I)) 
DIRNPCI+4)= DIRNS(I])+2.*DIRND(I)+2.*DIRNDD(T) 


TEAST. DIRECTION CCSING | 
DIRES(I)= DIREP (1) *ALPHAN«(DIREE(1)-DIREP(1)) 
DIRED(I)=DIRED(1-4)+2.*DIREDD(1-1)+BETAN/2,#(DIREE(I)-DIREP(I)) 
DIREDD(I)=DIREDD(1-1)+2.*GAMMAN/4.,«(DIREE(I)-DIREP(I)) 
DIREP(1+1)=DIRES(1)+2.*DIRED(1)+2,*DIREDD(1) 


DOWN DIRECTION COSINE 


DIRDS(I)=DIRDP(1I)+ALPHADe(DIRDEC(I)-DIRDP(I])) 
DIRDD(I)=DIRDD(1-4)+2.*DIRDDD(1-1)+BETAD/2,*(DIRDECI)-DIRDP(I)) 
DIRDDD(CI)=DIRDDD(I-1)+2.*GAMMAD/4,.*(DIRDE(IT)-DIRDP(I)) 

DIRDP (1 +4)=DIRDS(1)*2.*DIRDD(I)+2.*DIRDDD(I) 


COMPUTE SMOOTHED TARGET POSITION 


ZETACI)=ACOSF(DIRDS(CI))-PI/2. 
XTNMSCT)=XONMCT)+RNMSC IT) ®COSFCZETACI))*DIRNS(CI) 
YTNMSCI)=YONMCT)+*RNMSC IT) *®COSFCZETACI)) *®DIRES( I) 
ZTNMSCT)=ZONM(T)+DIRDSC I) *#RNMS(I) 


COMPUTE SMOOTHED AND ERROR DIFFERENCES FROM ACTUAL 


XTSMOCTI=XTNMCIT)-XTNMSCI) 
YTSMOCTI=AYTNMCT)<YTNMSCT) 
ZTSMOCT)=SZTNMCT)-ZTNMS(I) 
XTERRCI)=XTNMCI)-XTNMECI) 
YTERRCI)=YTNMCI)=-YTNME(I) 
ZTERR(I)=ZTNMC I) -ZTNMECI]) 
RSMO(1)=RNMCI)-RNMS(T) 


COMPUTE THE ACTUAL AND SMOOTHED SD RMS ERROR VALUES 


CIRERR(I)=SQRTFCXTERR( I) we2+YTERR( I) e*2+ZTERR(I] )**2) 
CIRSMOCI)=SQORTFCXTSMOCT) #e2+YTSMO(])*#*2+ZTSMO( I) **2) 


CHECK TO SEE IF THE ACOSF ARGUMENT HAS BEEN EXCEEDED, 
IF IT HAS SET THE COSINE EQUAL TO ONE AND PRINT OUT THE VALUE 


[IF (DIRNS(1)-41,.)2441,211,210 a 
PRINT 212,(T(1),DIRNS(I)) 

FORMAT (6X2HT= F8,4,2X6HDIRNS=F 8,4) 

DIRNS(1)=4.0 

[F (CDIRES(1)-4.)221,224,229 
PRINT222,(T(1),DIRES(1)) 

FORMAT (6X2HT=F8,4,2X6HDIRES=F8. 4) 

DIRES(1)=4.9 

IF (DIRDS41)744.0)251,2324,, tent o> yo gp logy a 
PRINT232, (T (bX, DIRDELTS (Pe is te 
FORMAT (6X2H7EF B44, SNGHDIRDSEFS » 4) ee 
DIRDS(1)=1,0 

CONTINUE 


CONVERT THE DIRECTION COSINES TO DEGREES FOR GRAPH OUTPUT 
enrMi Dias iat 
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120 


121 


25 
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ANSMOCT)=ACOSF (DIRN(CI))-ACOSFCDIRNS(I)) 
AESMOCI)=ACOSF (DIREC I) )-ACOSFCDIRES(I)) 


ADSMO(IT)=ACOSF(DIRD(I))-ACOSF(DIRDS(I)) 
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COMPUTE THE DESIRED LIMIT POINTS FOR THE GRAPH OUTPUTS 


(atl 


[FATECKS* WL) 2120 aeo eae? 
TCCK*#A STOW) FLU, 
AERRP(K)=ATANF(.5/RNM(1)) 
AERRM(K)'s=AERRP(K) 
RERRP(K)=,5 

RERRM(K)=-,5 
CERRP(K)=.875 

K=Ke4 

CONTINUE 


y @ 
i N Ew e ff ‘ ) ey ean 


IF DETECTION WAS MISSED INCREASE Q@ JN THE RDOT SMOOTHING 
COMPUTATIONS 


IF (ND(1))55,55,56 
Q=Q+1, 

GO TO 29 

=). 

TC] +4). =T(]) +2, 


END OF MAIN PROGRAM DO LOOP 
CONTINUE 
OUTPUT CONTROL 


COMPUTE MARKERS WHICH INDICATE THE TARGETS INITIAL POSITION 
ON THE SMOOTHED PROFILE CURVES 


X1(1)=XTONM+.2 
X1(2)=XTONM-,2 
YL(1L)=YTONM+.2 
Y1(2)=YTONM~,° 
Z1(1)=ZT0/6076,1+.09 
Z1(2)=ZT9/6076,1-.92 
X2(4)=XTONM+.2 
X2(2)=XTQONM-,.92 
Y2(1)=YTONM-,2 
Y2(2)=YTONM+,92 
Z2(4)=ZT0/6076.1-.09 
Z2(2)=27T9/6076,14+,.02 


GRAPH STATEMENTS 


IF (NGRAF1)279,279,271 
CONTINUE 

DO 191 I=1,4 

IT(1)=8H 

17(4)=8HGENTZ NE 
1T(2)=8H, PLANE A 


IT(S)=8HCTUAL:; +, 13 ; cli ; aah te 

1T7(4)=8HNCI SE. iC) mee) ah on der Dd dO 

5 

CALL DRAW(NUMPTS, YONM,XONM,1,0,9H 21T p42. 9120007 %r epee, One 
1LAST) 


CALL: DRAW(NUMPTS,YTNM,XTNM,2,054H 
CALL DRAW(NUMPTS, YTNME,XTNME,3S,0,4H 


) 
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Gos Ar4 
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www 


aaa 44 
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103 


104 


110 


270 
282 


105 


106 


116 


DO 103 1=1+4 TE 

1T(1)=8H atal 

IT(1)=8HGENTZ NE 

1T7(2)=8H PLANE §S 

1T(3)=8HMOOTHED 

174) =8HPOSITION 

lL =a > ee ie 

CALI. PRAWTNUMPTS, YONM. rou i, Tr 27T52,5,,19,,0,5,2,2,6,651, 
fi 





LLAST) 


CALL DRAW(NUMPTS,YTNMS,XTNMS»2+0,4H ) 
CALL DRAW (2,Y1,X4.2,0.4H ) 

CALL DRAW (2,Y2,X2036054H ) 

DO 104 1[=4,4 

IT(1)=8H 

1T(1)=8HGENTZ ND 

1T(2)=8H PLANE A 

IT(3)=8HCTUAL + 

17(4)=8HNOISE 

a 

CALL DRAWC(NUMPTS, XONM,ZONM,1,0,4H 91T049 09109200 50s Ge One 


1LAST) 


CALL DRAW(NUMPTS,XTNM»ZTNM,2,0,4H ) 
CALL DRAWC(NUMPTS,XTNME,ZTNME,35,054H ) 
DO 110 [1=1,4 

IT(1)=6H 

IT(1)S8HGENTZ ND 

LiC2) -Sh or SANE S 


1T($)=8HMOOTHED 
- eile ace 


ial DRAW(NUMPTS,XONM,Z0NM,1,0,4H olTst 9s oles 0s Ov arnesopo tues 


LLAST) 


CALL DRAW(CNUMPTS,XTNMS,ZTNMS,2,054H ~*~ ) 
CALL DRAW(2,X1,21,2,0,4H ) 

CALL DRAW (2,X2,22,3,054H ) 

CONTINUE 

IF (NGRAF2)281,281,282 

CONTINUE 

DO 4195 I1=4,4 

ITC 1)=8H 

ITCL)=BHGENTZ RA 

IT(2)=2=8HNGE SMOO - 
IT( 3) S8HTHING RE 

1T(4)=8HSULTS 


T=1 

CAGE DRAWCNUMPTS,7T,RSMO,1,054H LT eDdG code S00 0002 r6,8,1,LAST) 
CALL DRAWCNERMAX,TC,RERRP,2,2,54H ) n 

CALL DRAW(NERMAX,TC,RERRM,3,2,4H ) 

DO 106 [=1.,4 

IT(1)=8H 


IT(4)=S8HGENTZ NO 
IT(2)=8HRTH SMOO 
IT(S)S8HTHING RE 
IT(4)=8HSULTS 


IT=1 Totes er eas ' 

CALL LIRAV.INUMPTS3 T, ANSMO, 1,0; 44 , 59. Ot, 53,0,2,2,6,6,1,LAST) 
CALL’ DRAWINERMAAS TC>AERRPY20204H °°) ; 

CALL DRAW(NERMAX,TC,AERRM, 9,2, 4H ) 

DO 116 I=1.4 

1T(1)=8H 


IT(1)=8HGENTZ EA 
I1T(2)58HST SMOOT 
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IT(3)=8HHING RES 
IT(4)=8HULTS 


T=1 

CALL DRAW(NUMPTS,T,AESMO,1,0,4H e1Te 20. aeOleSe0 reser, Opie 
CALL DRAW(NEPMAX; [C,MERRE,2,254H 2) | 
Cat JRAW(NEPMAX, TC, AERRM,3,2,¢H | >" ss 


Dc eden ied 

108 IT(1)=8H | 
IT(1)=8HGENTZ DO 
1T(2)=8HWN SMOOT 
IT(3)=8HHING RES 
114) Gime 


CALL DRAW{NERMAX,TC,AERRP, 2, 2, 4H 
CALL DRAW{NERMAX,TC,AERRM,S,2,4H ) 
281 CONTINUE 
IF (NGRAFS)283,283,284 
284 CONTINUE 
Do 285 j=4,4 
285 IT(1)=8H 
IT(1)=6HGENTZ 3D 
1T(2)=8H RMS SMO 
IT(3)=8HOTHING R 
17(¢4)=8HESULTS 


bet . 
CALL DRAW(NUMPTS,T,CIRSMO,1,0,4H 21T, 500817200 + 0,6, ep 0,0me 
TEAST 

CALL DRAW(NERMAX,TC,CERRP, 3,2, 4H ) 


283 CONTINUE 
PRINT STATEMENTS 


PRINT 99 

95 FORMAT (1H1,27X17HSMOOTHING RESULTS) 
PRINT 64, °¢17(1),125,6) 

64 FORMAT(/s,352X,2A8) 
PRINT 107 


CALA DRAW(NUMPTS,T,ADSMO,4,0,4H 1 1T 15004 0415000202060 6e1e LAST) 
) 


107 FORMAT( /,6X66HT2TIMEs R NMS=ACTUAL RANGE, CIRERR=SDIMENSIONAL RMS 


4 ACTUAL ERROR, »/,OX64HCIRSMO=3 DIMENSIONAL RMS SMOOTHED ERROR 
2, ND=1 DETECTION, ND=0 NO,/,6X 68HDETECTION, NG=QUALITY COUNTER, W 
mriel NQ=9 TARGET CONSIDERED LOST AND ,/7,6X24HALL NINES ARE PRINTED. 
) 
PRINT 96 
96 FORMAT ( /st5X4HT,4X4HR NM, 7X6HCIRERR,4XOHCIRSMO,2X2HND, 
12X2HNQ) 
PRINT 97,¢T(1)/RNM(I)-CIRERR(I)/CIRSMO(T),ND(1)sNQ(1),1=1,NUMPTS) 
97 FORMAT ( 7X1F10.0,35F10.4,214) 
2 CONTINUE 
400 STOP 
END 
END 
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Appendix’ C 


Computer Program of Program ALPHA CHECK 


Only the smoothing and prediction and output sections are shown 


since the preceeding portions are the same as in THREE D. 


, ’ ’ t ’ t ¢ 4 
; a / vei Ny * ’ ‘ & 
Pe es Aig = o7]l of Loo ' Hoa 
' , ‘ 1 1 ww ' 
car ? 4 4 ’ t 


{i Ta i an : i e f v »Pe 





47 IF (RNM(1)-80.)65,66,66 
65 PD=1,55-.016*RNM(]) 
GO TO 52 
66 IF (RNM(I)-110,)62,68,68 
62 PD=.88-.0077*RNM(J) 
GO%t Oe 22 4 toe 
68 FD=.297,001¢*RNA()) >. . ) o) rary a 4 
52 CONTINUE 
ADEV=ABQSF (DEV) 
ADEV=0.° 
IF (PD-ADEV)53,54,54 
53 NOCI)=0 
GO TO 29 
24 NQ(I)21 
20 TCI+1)=TC1) +2, 
1 CONTINUE 





SMOOTHING COMPUTATIONS 


FPTS=NUMPTS 
NALPH=20 


RANGE SMOOTHING CONSTANT ANALYSIS 


ALPHAR=,95 
DO 304 J=1,NALPH : 
P=4. 
AR(J)=ALPHAR 
ERAR(J)=Q. 
DO $95 [=1,NUMPTS 
IF (NQ(1))502,592, 593 

302 IF(P)304,304,305 

394 RNME(I)=RNMP(T) 

7 RDOTME(I)=RDOTMP(]) 
GO To 306 

SoS I1F(P)3596,396,598 

308 P=P-1. 
eie=¢1)/2.+1. 

RNMP( I) SRNMECT) 
RDOTMP(I)=RDOTME(1) 
RDDOTE(I-1)=0. 

306 RNDOTE(])=RDDOTE(1-1)+1.7(2.4%Q)*(RDOTME(])*RDOTMP(])) - 
RNMS(1)=RNMP(1)+ALPHAR*(RNMECI)-RNMP(I)) 
RNMP(1#41)=RNMS([) +2. #RDOTME(])+2, *RDDOTE(!) 
RDOTMP(1+4)=RDOTME(I)+2.*RDDOTE(1) 
RSMO(I)=¢(RNM(1)-RNMS(I)) «%2 
ERAR(J)=ERAR(J)*RSMOC(CI)/(FPTS“CTC) 
IFCNQ(])) 55,55,56 

55 Q=Q+1. 
GO To 305 
SOMN=1. 
$05 CONTINUE 
301 ALPHAR=ALPHAR+,05 


T=1 
PRINT 414 >)» 28°: ) 2 ) Fos ,OO ) ) 
114 FCRMAT(GH1G47¥O6KOPTIMUM, SHOCGTHING CONSTANT ANALYSIS,//,6X62HQUAL 
1TY INDICATOR, THE MINIMUM AVERAGE MEAN SQUARE ERROR ) 


PRINT 604,°0IT(I),1=5,6),ANOISE) 

604 FORMAT(//,20X.2A8,2X7HANOISE=F 49.2) 
PRINT 115 

115 FORMAT(//,29X16HRANGE, ALPHAR ) 
PRINT 116, (AR(J),ERAR(J),J2=1,NALPH) 


? OO Se 





116 FORMAT(//,((23X7HALPHAR= F4.2-2X5SHERAR= 


322 
524 


323 
528 


$26 


$25 
s24 


iy 


0 


552 
554 


5338 
558 


556 


535 
551 


ArrULNWILA WY 


F6.4)/)) Eien ai 


NORTH DIRECTION COSINE SMOOTHING CONSTANT ANALYSIS 


ALPHAN=.95 

DO 223i K=4/NALPH 
Po et 
ANC K)SALPHAN 
ERAN(K)=0, 
BETANZSALPHANww2/(2, -ALPHAN) 
GAMMAN=2. *BETAN/ALPHAN®@=(ALPHAN*+BETAN 
DO $25 JT=4,NUMPTS 
IF(NQ(J))$22,322,323 

IF (P)324,324,325 

DIRNE(J])=DIRNP(I) 

GO TO 326 

IF (P)326,326, 528 

P=Pe-1. 

DIRNP(I)=DIRNE(]) 

DIRND(I=1)=0. 

DIRNDD(I-41)20. 


a 


§ ’ out #4 , + b 


) 


DIRND(T)=DIRND(1-4)+2.*DIRNDD(1-1)+BETAN/2,*(DIRNE(I)-DIRNP(I)) 


DIRNDD( JT) =DIRNDD(I-1)+2. *GAMMAN/4.*(DIRNE(I)-DIRNP(I1)) 
DIRNSCI)=DIRNP( 1) +ALPHAN*(DIRNECI)-DIRNP(1)) 


DIRNP(CT+*41)=DIRNS(I)+2.*DIRND(I)+2.*DIRNDDCI) 


DNSMO(1)=(DIRN(I1)=DIRNS(1)) «#2 
ERAN(K)=ERAN(K) #DNSMO(])/(FPTS-CTC) 
CONTINUE 

ALPHAN=ALPHAN+,05 

[=4 

PRINT 147 


PRINT 118, (AN(K),ERAN(K) ,K=1,NALPH) 


FORMAT(//,20X34HNORTH DIRECTION COSINE, ALPHAN ) 


FORMAT (//s ((2SX7HALPHAN= F4.2,2X5HERAN= F10.8)/)) 


EAST DIRECTION COSINE SMOOTHING CONSTANT ANALYSIS 


ALPHAE=,.05 

DO 534 L=4,NALPH 

P=i. 

AE(L)=ALPHAE 
BETAE=ALPHAEx#2/(2,-ALPHAE) 
GAMMAE=2. *BETAE/ALPHAE~(ALPHAE+BETAE) 
ERAE(L)=0, 

DO $35 J[=4,NUMPTS 

IF (NQ(1))332,332,333 

IF (P)334,334,335 
DIREP(I)=DIREEC(T) 

GO TO 336 

1F(P)$36,3356, 5358 

P=P-1, 

DIREP(I)=DIREE(!) 
DIREDCi<-2)=07 
DIREDD(I-1)=0. 


DIREDK i s2DIRED( F-41542. eDiREDD( 1-1: +8CTAEs/2.x(DIREE(I)-DIREP(I)). 
DIREDI (1) =DIREDE (1-142. *GAMMAES4 #(DIFER(*) @-DIREP(I)) 
DIRESCLIEDI REP (1) +aALPHAEw(DIREE( Ls -G1REP( I>) 
DIREP(CI#1)=DIRES(1)+2.*DIRED(1)4+2,*eDIREDD( 1) 


DESMOCT)=(DIRECT)-DIRES(I))**2 
ERAE(L)==ERAE(L)+DESMO(I)/(FRPTS-CTC) 
CONTINUE 

ALPHAE=ALPHAE+,05 


«- he “ “ ee al 








119 


161 


342 
344 


3438 
548 


546 


545 
541 


175 
123 
274 


110 
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fe 
PRINT 119 
FORMAT(1H1,20XS1H EAST DIRECTION COSINE, ALPHAE 
PRINT 161, (AE(L),ERAE(L)»L=1,NALPH) 
FORMAT(//, ((23X7HALPHAE= F4.2,2X5HERAE= F10.8)/)) 

; DOWN, DIRECTION COSINE SMOCTHING GONSTANT ‘ANALYSIS 


>>> Waa POU NN Ne wee sees 


ALPHAD=,.05 
DO 341 M=4,NALPH 

P= ie, 

AD(M)=ALPHAD 

BETAD=ALPHAD«#2/(2,-ALPHAD) 

GAMMAD=2. *¥BETAD/ALPHAD=(ALPHAD+BETAD ) 

ERAD(M)=Q, 

DO 345 J=4,NUMPTS 

IF(NQC 1) 342,342,343 

1F(P)344,344,345 

DIRDE(I)=DIRDP(]) 

GO TO 346 

IF (P)346,346,348 

P=P-4, 

DIRDP(I)=DIRDE(]) 

DIRDD(I-1)=0, 

DIRDDD(I-41)=9, 

DIRDD(I)=DIRDD(1-1)+2.*DIRDDD( 1-1) +BETAD/2. *(DIRDE(1)=DIRDP(I)) 
DIRDDD(CI)=DIRDDD( I -1)+2. *GAMMAD/4. #(DIRDE(1)-DIRDP(1)) 
DIRDS(T)=DIRDPC(I)+ALPHAD*(DIRDE(1)-DIRDP(I1)) 
DIRDP(I+1)=DIRDS(1)+2.*DIRDD(I)+2.*eDIRDDD(1) 
DDSMO(I)=(DIRD(1)-DIRDS(1)) «2 
ERAD(M)=ERAD(M)+DDSMO(I)/(FPTS-CTC) 

CONTINUE 

ALPHAD=ALPHAD+,05 

l=1 


mer INt 122 


FORMAT(//,20X3S1LHDOWN DIRECTION COSINE, ALPHAD  ) 
PRINT 123,(AD(M),ERAD(M),M=1,NALPH) 
FORMAT(//,((23X7HALPHAD= F4.2,2X5HERAD= F10.8)/)) 
IF (NGRAPH)279,270,271 
CONTINUE 
DO ede t) l21,4 
IT¢(1)=8H p 


IT(LI=8HGENTZ AL 

IT(2)=S8HPHAR VS 

1T(5)=8HAVG MEAN 

IT(4)=8HSQ ERROR 

CALL DRAW(NALPH,AR,ERAR,0,0>4H 211s M0. + Or0% eres 5,6, 1, 0KGTe 
DO 111 1=4,4 

IT(1)=68H 

IT(CL)=8HGENTZ AL 

IT(2)=8HPHAN VS 

IT(3)=8HAVG MEAN 

1T(4)=8HSQ ERROR 

es. ©, «oe ey hs a 
CALL DRAWOCNALPR YAN, ERAN, O,90,44 
DO 112 [=4,74> : : | | 
IT¢(1)=8H 

IT(1)=8HGENTZ Al 

TT(2)=8HPHAE VS 

IT(S)=8HAVG MEAN 

1T(4)s8HSQ ERROR 


) 


91744256 60050%252,5,6,1, LAST) 


qa, 2 & see ee eee _eeee aoe lCU CCl 





CONFIDENT RA® 


113 


270 


APPENDIX C 
La 
CALL DRAW(NALPH,AE,ERAE.Q,04H 1 Tr elo ee Or Qed nes cisGpiie none 
BO tio Tat, 4 
[T41)=8H, She ae i a ae 
DTXet ) SBC ANT 2 AL Se ee 
ITCA)SB8HPHAD VSSQ totes tee eee eee eee vee 
1T(3)=8HAVG MEAN 
17(4)58HSQ ERROR 
lead 
CALL DRAW(NALPH,AD,ERAD,0,0,.4H r1To e200 00 rQeeoer2, Goto bAST) 
CONTINUE | 
ANOISE=ANOJSE+¢,5 
CONTINUE 
END 
END 
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SMOOTHING RESULTS 





RUN 2A 


T=TIME, R NM=ACTUAL RANGE, CIRERR=3DIMENSIONAL RMS ACTUAL ERROR, 
CIRSMO=S -DIMENSTONAL RMS SMOOTHED. ERRORy NDA DETECTION, ND20 NO 
DETECTION,: NI=QUALITY: COUNTER,’ WHEN: NQ=0; TARGET CONSIDERED LOST AND 
ALL NING@S ARE PRINTED, » bs 25> 1h > 


7 9» ope > o>» 


il R NM CIRERR CIRSMO ND NQ 
14,1508 ,2429 2929 q 2 
44,5361 1,8578 1884 0 1 
44,9288 510074 2,53538 1 5 

~—©45,3280 2,6577 2.4972 it 5 Ui 
19 7oa7 2591 1.9092 4 7 
16,1559 5,243 9 2D 6 ee 
6.9840 ,2384 41.3833 1 10 
17,9281 2,8164 , 3864 1 10 
17,4867 , 2849 »4138 1 10 
17,9446 3,1760 1.2653 1 
18,4455 . 9857 ~6776 ye ale 
18,8975 14,9631 gO 7 wae mag 
19,3867 ,2697 29958 4 140 
19,8789 5462 ~2026 1 10 
29.3618 1.0472 649439 4 49 
20.8912 4633 ., 49537 Jel 
21,9454 , 8665 .6597 4 10 
24,8344 $.2078 4.2478 0 
22,3537 2.35855 359 57 1 10 
22.8417 ~6766 2586 et 
25,3525 2.2829 1.97356 1 10 
25,8669 ~4162 7837 oo 
24,3896 8367 4905 4 10 
24,9098 2,7580 1.18538 1 10 
48, 25,4187 14,8427 .4945 14 10 
Bis eo 2U7 14,8184 68744 1 10 
52, 26.4281 2.0666 1.9061 1 10 
54, 26,9446 .1707 .7198 aL ole 
D0. 27,4663 ~9916 9665 LP 10 
58, 27,9841 2.4726 14,2136 4&4 10 
60, 28,4961 22,8985 1.9486 eo : 

62, 29,0163 ,9102 68421 1 10 
64, 29,5499 2.4250 15D S45 1 10 
66, 30,0547 .2588 7/57 1 10 
68, 39,3636 2.8572 4.9178 4 10 
70. $1,07352 41.1682 7910 Ta 4 0. 
72.6 w1,5752 22,9874 ©9229 1 10 
74, $2,0719 4.1586 9257 0 9 
76, $2,5737 4.3505 ©9190 0 8 
78, 53.0744 73983 - .29350 1 10 
80, 53.5744 4837 ~2606 0 
82, 54,0735 3.1434 1.3395 ie 0 
84, 34.5622 6.2619 2.0776 1 10 

86,'" SawURS7 =.) S705,” eos oe ae 

BB, 2 Oo aS Soon S eco en 

90,  o9-9716 » 2 oe2o1 aeenieriey ot sated) a 
92. 36.4471 ~8415 9765 1 10 
94, 56.9295 4.1639 10170 0 9 
96, 57,4029 .7998 1.0545 0 8 
98, 37,8898 2.9077 4.8114 14 10 
100. 58.3598 2.2822 1.9747 oe 
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